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(57) Abstract 



A method is provided for rapidly producing high resolution two-dimensional and three-dimensional real-time images of an object 
by MRI imaging comprising, acquiring a set of data for the object at a raster of coordinate points in a region in a k-space using a first MRI 
imaging technique and updating some of the data using a second MRI imaging technique. 
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KEYHOLE MRI 
FIELD OF THE INVENTION 

The invention relates to magnetic resonance imaging techniques and in particular to 
magnetic resonance imaging techniques for rapidly producing high resolution two dimensional 
5 and three dimensional real time images. 

BACKGROUND OF THE INVENTION 
Magnetic resonance imaging, "MRI" is an important noninvasive methods for looking 
into the human body. Medical MRI scanners first became available in the early 1980s and since 
then MRI equipment and methodology has advanced rapidly. Today, conventional MRI 
10 scanners are capable of imaging internal body structures and resolving details smaller than a 
cubic millimeter in size and MRI scanners with microscopic resolution a thousand times better 
(0.1 mm cube ) are available. 

MRI is based on the phenomenon that when a sample material, including human tissue, 
is introduced into a magnetic field B, nuclei in the sample having intrinsic spin will precess 
15 around the direction of B with a component of spin parallel or anti-parallel to B. This is 
because nuclei with spin behave in many ways like small dipole magnets that spin around their 
own axes. The precession frequency about the magnetic field is known as the Larmor 
frequency. The relationship between the nuclear spin, 5, and the magnetic moment, //, of a 
nucleus is /}= y* S where the proportionality constant, y, is the gyromagnetic ratio of the 
20 nucleus, which is distinctive for each nucleus. The Larmor frequency, od, with which the spin 
and associated magnetic moment of the nucleus precess around the magnetic field is equal to 

yB. 

Slightly more of the nuclei are aligned with components of spin parallel to the direction 
of B than are aligned anti-parallel to B because nuclei aligned with B are in a lower energy 
25 state than nuclei aligned against B. The material is therefore polarized with a net spin density 
per unit volume, f{x,y,z), and a magnetization density per unit volume, m(x,y,z) = 
y* jo(x, y, z) , pointing in the direction of B. 

The gyromagnetic ratios for some common nuclei are.-lR, 4257 Hz/gauss; 13c, 1070 
Hz/gauss; 23Na, 3368. In a 1.5 Tesla polarizing B field these gyromagnetic ratios resuh in 
30 precession frequencies of 63.85 MHz, 16.05MHz, and 50.5 IMHz respectively. For a given 
polarizing B field, however, the obser\'ed Larmor frequency for a particular nucleus is shifted 
from the expected value yB and depends upon the molecule in which the nucleus is found and 

1 
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how the nucleus is bonded in the molecule. The shift is caused by currents that are induced in 
the electron cloud of the molecule by B that shield the nucleus from B. The magnitude of the 
shift is substantially proportional to B and can be used to identify the chemical environment in 
which the nucleus is found. This shift is called a chemical shift. 
5 The Larmor frequencies of nuclei, because they can be used to identify the chemical 

environment in which nuclei are found and because they are dependent on the polarizing 
magnetic field B, can be used to effectively map the density distributions of nuclei and 
molecules in an object. If the densities of the nuclei and molecules vary with different internal 
structures of the object then these internal structures can be identified and imaged. 

10 The basic approach to using Larmor frequencies for imaging an object is to produce a 

magnetic field that changes over the volume of the object so that the Larmor frequencies of 
nuclei investigated var\^ as a fimction of position. The nuclei are then stimulated to radiate 
energy at their Larmor frequencies. The signal intensity at a particular frequency is 
proportional to the number of radiating nuclei at the position in the object which is identified 

15 by the signal frequency. As a result, mappings of the densities of various nuclei in the object 
can be made which are functions of position in the object and intemal structures of the object 
can be visualized. 

Implementing the above procedure is generally performed by polarizing the object with 
a strong homogeneous magnetic field, Bq, in the z direction. Bq polarizes (along the z axis), all 

20 nuclei having spin and magnetic moment, and each type of nucleus precesses around Bq with 
its own Larmor frequency. An imaging nucleus, i.e, nucleus whose density it is desired to 
image having a gyromagnetic ratio y and Larmor frequency co 0=7^0' selected for imaging 
by rotating its spin density p(x,y,z) and associated magnetization density m(x,\\z) away from 
the z axis. This produces a net spin density and magnetization density perpendicular to z which 

25 rotate with Larmor frequency coq- The rotation is achieved by applying an RF selection pulse, 
having a narrow frequency spectrum centered at the Larmor frequency coq of the imaging 
nucleus and magnetic field Bgx along the x axis, for a short period of time Tg. Bgx causes the 
net spin density of the imaging nucleus to rotate away from the z axis by an angle yBs^Ts^ 
conventionally called the flip angle. A receiver coil positioned at the side of the object with the 

30 plane of the coil parallel to the z axis (z.e. the perpendicular to the plane of the coil is 
perpendicular to the z axis) senses the rotating magnetic field moving with the rotating 
magnetization density. The rotating magnetic field illuminates the receiver ce?Q times a second, 
generating a signal, SIG(t), in the coil, which oscillates with frequency coq. 
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Gradient fields in the z direction, which vary hnearly with x, y and z, are then appUed to 
the object. The gradient fields cause the frequency of rotation of the magnetization density to 
vary around Larmor frequency coQasa function of position in the object. If Gx, Gy, and G^ are 
the derivatives of the gradient fields, the signal becomes, to within a constant: 

5 SIG(t)= JJJ/7(x, \-,z)exp[zY(Bo+GxX+Gyy+G2z)t] JxiivJz= exp[zcoot]S(t), S(t) is the position 

dependent part of the integral, it has been assumed that the angle of rotation of p{x,y,z) away 
from the z axis is 7r/2, and relaxation of the spin density to zero in the transverse plane and to 
equilibrium along the z axis is ignored. 

By changing variables so that kx= (YGxX)t, ky= (yGyy)t and k^^ (7G2z)t, S(t) can be 

10 written as a function of position in a ''k-space". S(t)^ S(kx,ky,k2)^ S(k), and: 

S( k )= JjJ/?(^% y'y^) exp[z (kxX+kyy+kzZ)] dxdydz = \\l p{r )exp[i k • r jd^r. 
The last integral is the Fourier transform of the spin density function p{x^y\z) of the 
imaging nucleus. S(k) is a *'k-transform'' of p{x,y,z) and k and r are conjugate variables so 
that /7(r)-JJJS(k)exp[ik-r]d3k. 

15 When the imaging nuclei are first flipped away from the z axis they precess together 

coherently with a net spin density and magnetization density in the xy plane. With time, 
however, the coherence in the transverse xy plane decays to zero and the spin density relaxes to 
the equilibrium state where the imaging nuclei are polarized along the z axis. The decay of net 
transverse spin density and return to equilibrium along the z axis are characterized by different 

20 time constants known as T2 and Tj respectively. When inhomogenities in the magnetic field 
are present the decay of transverse spin density is accelerated and is characterized by a time 
constant known as T2*. The relaxation times are related by the inequality T2* < T2 <Ti . 

The various relaxation times of imaging nuclei are quite sensitive to the chemical 
environment of the imaging nuclei and the molecular structure of the molecules in which they 

25 are bound. This sensitivity allows different materials to be discriminated from each other by the 
rate at which their MRI signals decay. In particular in medical and biological imaging different 
tissues and features in biological tissues can be identified and contrast enhanced by 
manipulating their T^, T2 and T2* time constants. T2* for example is of particular interest in 
functional neuroimaging. 

30 Often the three dimensional imaging problem is reduced to a two dimensional one by 

selecting nuclei for imaging that are confined to a thin transverse "xy" slice of an object 
imaged. Assume that while the RF selection pulse is on, a selection gradient field that varies 
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linearly with z with a gradient Ggz is applied to the object parallel to the z axis. The Larmor 
frequency of the imaging nuclei selected will then be shifted ft-om cvq as a, function of position 
along the z axis by an amount ^ Gg^z . If the RF selection pulse is a narrow pulse with a 
central frequency co g, only those imaging nuclei in a thin transverse slice at z coordinate Zg 
5 such that C(?^= cvq will be rotated away from the z axis by the RF selection pulse. 

By acquiring data for two dimensional k-functions, S(kx,ky,Zs), for a range of values Zg, a 
three dimensional image of the object can be constructed. 

Many different techniques and procedures have been developed for MRI imaging. All 
involve procedures for evaluating the k-transform, S(k), of an object at many points in a raster 
10 of points in k-space so that the Fourier transform of S(k) results in a proper evaluation of 

Typical MRI imaging techniques are the conventional gradient echo (GE) and spin echo 
(SE) imaging sequences. Each of these sequences acquire k-transform data for a slice of an 
object imaged at points in k-space (hereafter "read points") along lines in k-space (hereafter 

15 "read lines") parallel to the k^ axis. The sequences begin with an RF slice selection pulse and 
acquire data for only a single read line each time they are run. To acquire data for a complete 
slice of the object the sequences are repeated as many times as necessary to read data along all 
the read lines required for the k-transform of the slice. The sequences are similar except that the 
SE sequence includes a focusing or rephasing RF pulse for neutralizing phase differences in the 

20 spins of imaging nuclei arising from magnetic field inhomogenities. 

GE and SE sequences can be used to acquire high resolution k-transforms using 
moderate gradient fields and receiver bandwidths. However, acquiring data using these 
sequences is relatively slow. Delay times between repetitions, Tr, of both these sequences are 
generally long. Tr depends upon, among other factors, how fast the spin density which is 

25 rotated away from the polarizing field recovers its equilibrium value along the polarizing field 
direction after a repetition of the sequence. Tr must therefore be longer than the time, Tj, 
which is characteristic of the recovery if repeated sequences are all to have the substantially the 
same signal strength. 

In most plant and animal tissue, T^ ranges from 0.2 to 1 sec. Therefore, conventional 
30 GE and SE sequences are inherently slow, with total image times for a slice ranging from 
seconds to minutes. However, their long imaging times restrict their use in real time rapid 
imaging of motion of internal organs and metabolic changes in the human body. 
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Real time imaging of physiological processes such as cardiac contractility, peristalsis 
and brain activity are important areas of medical research. MRI techniques are needed for high 
resolution rapid imaging for these and other dynamic biological processes. However, even the 
most basic of the body's internal movements and metabolic changes requires imaging 
5 techniques that can produce clear accurate pictures in tenths of a second. For example, the heart 
beats with a frequency on the order of one beat per second. To follow the progress of a single 
heartbeat, or assess an arrhythmic event, it would be advantageous to acquire high definition 
images of the heart on the order of every 0.1 seconds or less. 

Techniques exist for fast temporal MRI imaging. One such technique is echo planar 
10 imaging (EPI). In EPI, all or a large part of the area of k-space of a slice to be imaged is 
scanned in a short period of time following a single RF slice selection pulse. Data for imaging 
a complete slice of an object using EPI can be accumulated in a time period of from 20 to 100 
milliseconds. While these time frames are useable for imaging physiological processes, the 
images obtained with EPI are generally inferior in resolution to images produced using 
15 conventional GE and SE techniques. 

Fast EPI techniques also require special hardware and large receiver band widths and 
they use very strong rapidly changing gradient fields. The rapidly changing fields can produce 
undesirable biological effects, such as electrical stimulation of nerves and slight torques on 
cones in the retina which are evidenced by light flashes in the eye. 
20 Fig. 1 A shows an EPI sequence 20, for acquiring with one repetition of EPI sequence 

20, a two dimensional k-transform, S(kx,ky,Zs), (hereafter, S(kx,ky) ) for a complete shce of 
an object at N x M read points. An RF shce selection pulse 24 and a selection gradient field 22 
and 26 are applied to the body being imaged to select a slice to be imaged at a z coordinate Zg. 
Data is taken using a phase encoding gradient 28, comprising N-1 phase blips 36, and a read 
25 gradient 32, comprising N read blips 38, of constant magnitude and alternating sign. Read blips 
38 are Tjyj seconds long. 

Phase encoding gradient field 28 begins with a negative gradient 30, which dephases 
the spins and sets ky to a minimum negative value. This minimum negative value is the ky 
coordinate of a first read line of M points in k-space at which data is sampled. A read gradient 
30 field 32, begins with a negative gradient 34, which dephases the spins and sets kx to a 
minimum negative value which is the kx coordinate of the left most point of each read line in 
k-space. Following gradients 30 and 34, read gradient 32 progresses with N read gradient blips 
38 during which data for S(k) is recorded . At the end of each read gradient blip 38 a phase 
gradient bhp 36 increases ky by a constant amount. Phase blips 36 and read blips 38 are 
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carefully timed so that the centers of phase blips 36 occur at the zero crossover points between 
positive and negative read blips 38. This timing ensures that each read line in k-space begins 
and ends with the same coordinates as the other read lines in k-space. 

A data signal 40, comprises a pulse train of N data echo pulses 42, one for each read 
5 blip 38. Each data echo pulse 42 is substantially centered on a read blip 38, and is sampled M 
times, usually at equal time intervals, during the duration Ty[ of each read blip 38. Therefore, 
with each read blip 38, a set of M data values for S(k) is acquired at M equally spaced points 
on a read line in k-space, parallel to the kx axis. EPI sequence 20 therefore records in one 
repetition N x M data values for S( k) in a total imaging time for the slice Ti= N x T^ seconds. 

10 A sampUng trajectory 50 in k-space, traversed by using EPI sequence 20 for sampling 

data at N X M points 51, is shown in Fig. IB. Trajectory 50 is a series of equally spaced read 
lines 52, parallel to the kx axis connected by short phase lines 54, parallel to the ky axis. The 
length of phase lines 54 is equal to the time integral of a phase blip 36, multiplied by the 
gyromagnetic ratio of an imaging nucleus whose density distribution is imaged. The length of 

15 read lines 52 is equal to the time integral of a read blip 38, multiplied by the gyromagnetic ratio 
of the imaging nucleus. Phase gradient 30 and read gradient 34, displace the beginning of 
trajectory 50 to point 56. Arrows 58, indicate the direction of traversal along trajectory 50 as 
data is sequentially taken at points 51 on trajectory 50. 

For rapid, snapshot imaging of motions and processes in the human body, conventional 

20 MRI imaging techniques have generally satisfactory resolution but are too slow. Fast 
techniques, such as EPI, are rapid enough, but lack resolution and have serious technical risks 
and drawbacks. 

It would be desirable to have an imaging technique that would have the resolution of 
conventional GE and SE techniques but produce images in time frames approaching that of fast 
25 EPI techniques. 

It is further desirable that such imaging techniques not require pulsed gradient fields 
with time rates of change that might produce undesirable biological effects. 

SUMMARY OF THE INVENTION 
It is an object of one aspect of the present invention to provide an MRI method for two 
30 and three dimensional MRI imaging that can be used to produce a sequential series of high 
resolution images of an object in periods of time short enough so that the method can be used 
for rapid, real time imaging, of processes in the body and motion of body organs. 
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A preferred embodiment of the method comprises acquiring a high resolution primary 
MRI k-transform, So(k), of an object at time tQ and updating only some of its values while 
retaining the rest unchanged to produce sequential k-transforms of the object, Si(k), S2(k),... 
Sn(k), at subsequent times ti,t2 . . .tj^, as the object changes. By updating only some of the data 
5 for each sequential k-transform it is possible to produce the sequential k-transforms very 
rapidly. By using data from the high resolution primary k-transform in the sequential k- 
transforms, it is possible to substantially preserve the high resolution of the primary k- 
transform in the sequential k-transforms. 

The values of So(k) which are updated are, preferably, only those that change 
10 significantly with the changes in the object. Those values of So(k) that do not change 
significantly with the changes in the object do not have to be updated and therefore are 
preferably retained from k-transform to k-transform. As a result, each of the k-transforms 
Sn(k), at times subsequent to the acquisition of the primary k-transform, has substantially the 
same number of accurate data points as the primary k-transform So(k). Therefore, the real 
15 space images of the object which are constructed from Si(k), S2(k),... Sn(k) at times ti,t2 
...tji all have approximately the same spatial resolution as the primary real space image of the 
object. 

The possibility of advantageously differentiating between values of So( k ) that change 
and values that substantially do not change with changes in an object imaged, is a result of the 

20 form of the dependence of So(k) on k. So(k) generally changes significantly only for points 
k in a limited volume of k-space of size l/L^, centered at the k-space origin, where is a real 
space volume characteristic of features of interest under observation in the object imaged. For 
points k outside k-space volume l/L^, So(k) does not change significantly. 

This dependence of So(k) on k can be understood from the relationship between the 

25 conjugate variables k and r . Let Rx be the x axis spatial resolution of an image constructed 
from a k-transform S(k) of the object. Let be the range along the kx axis over which data 
for S(k) is accumulated. Then, from a well known characteristic of Fourier transforms, 
KxRx= 1. Therefore, if an MRI image of an object is required to have a resolution R, in all 
three dimensions (for simplicity it is assimied the resolution is the same along all three axes 

30 though there is no necessity for this to be so), the primary k-transform of the object, So(k), is 
evaluated over a primar\^ raster of points in k-space for which -1/(2R) < ] k | < 1/(2R). 
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From similar considerations, objects, or features of an object, which have characteristic 
length of magnitude L, have significant values for So(k) at points k for which, 
|k| < 1/(2L). For changes on the order of ±AL, large relative changes in the amplitude and 
phase of So(k) occur in a bandwidth of k values for which i k | is in the region ko(l±AL /L), 
5 where kQ =1/L, and it is assumed that L is significantly larger than AL. For |k|»l/[2(L-AL)], 
essentially only the phase of So(k) changes while the magnitude remains unchanged. A change 
in phase is equivalent to a displacement in real space. Neglecting these phase changes for large 
values of I k I ~ 1/R results in a potential error of ~ 1/2R in the location of details of size 1/R in 
the real image of the object reconstructed from SqC k ). 
10 As a result, if changes of order ±AL in an object or feature of characteristic length L are 

being observed over the course of a time study, it is generally sufficient to update SqC k ), only 
for values of k for which | k | < 1/[2(L-AL)] in order to follow the changes with spatial images 
having substantially the resolution R of the primary spatial image. 

Therefore, in accordance with a preferred embodiment of the present invention, updates 
15 of So(k) at times ti,t2 . . .t^ are preferably taken at a raster of points bounded in k-space by a 
cube having edges of length on the order of 1/L which is centered at the origin of k-space. 

Values of So(k) at points in k-space outside this cube, are preferably retained from 
k-transform to k-transform. 

The number of points at which updates of the primary k-transform Sq( k ) are made, in 
20 comparison to the number of data points in the primary image, can be approximated by the 
ratio of the k-space volume of the update raster of points to the k-space volume of the primary 
raster of points. This ratio is equal to R^/L^ and r3/l3 for two and three dimensional k- 
transforms respectively. It can be seen that even small ratios of R to L lead to very large 
reductions in the number of data points in an update compared to the number of data points in 
25 the primary k-transform So( k ). 

Quite generally, for a feature in an object to be observed adequately, the resolution must 
be smaller than the feature and preferably significantly smaller than the feature (Le, R«L). 
The number of values of So(k) which change significantly with changes in the features of the 
object under observation is therefore often much smaller than the total number of data points in 
30 So(k). Updating can therefore generally be done very quickly in comparison to the acquisition 
time of the original primary k-transform So( k ). 
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This is not only because there are simply fewer data points to be evaluated for an update 
than for a primary image. The main advantage is that the reduced number of data points in the 
update permit the use of very fast MRI imaging techniques, such as for example EPI, to 
perform the updates. These techniques are not particularly attractive for rapid MRI imaging 
5 when used to acquire large amounts of data for a high resolution k-transform. When used to 
acquire data over large volumes of k-space they are expensive, require special equipment and 
very high, potentially harmful, gradient field time derivatives, and produce results with 
generally inadequate resolution. But they can be used to advantage for acquiring update data 
over a reduced number of data pomts. With the reduced number of data points these techniques 

10 can perform k-transform updates in time spans of from 20- 120msec with the moderate field 
strengths and receiver bandwidths used in conventional GE and SE techniques. 

It is thus seen that by sectioning k-space, in accordance with a preferred embodiment of 
the present invention, it is possible to use combinations of high resolution (for example, 
preferably, GE and SE sequences) and very fast MRI imaging techniques (for example, 

15 preferably, EPI sequences) to provide real time imaging of an object that is both very fast and 
very accurate. Furthermore this can be done without recourse to problematic gradient field 
derivatives or special equipment. 

It is realized that different MRI imaging sequences, and not only those mentioned 
herein, can advantageously be used for acquiring data over k-space regions that are sectioned in 

20 accordance with a preferred embodiment of this invention. It is therefore possible to combine 
different MRI sequences in accordance with a preferred embodiment of the invention to 
provide MRI imaging techniques that are sensitive to different types of temporal changes in the 
body, such as motional changes due to peristalsis or cardiac contractility or physiological 
changes resulting from neural activity. For example in some preferred embodiments of the 

25 present invention, preferably, GE or SE MRI imaging sequences are used to acquire a high 
resolution primary k-transform of an object imaged and updates are performed using, 
preferably, a fast EPI sequence. Since EPI is inherently T2* sensitive, MRI methods in 
accordance with a preferred embodiment of the present invention, are particularly suitable for 
use with contrast agents such as Gd-DTPA in perfusion studies, and in functional 

30 neuroimaging using contrast agents or dependent upon changes in blood oxygenation. By 
incorporating 180 degree RF pulses in the MRI sequences the methods can be made T2 
sensitive instead of T2* sensitive. 

In other preferred embodiments of the present invention an MRI imaging technique 
comprises a diffiision sequence so that a diffusion weighted MRI image of a subject can be 

9 
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obtained. Diffusion weighted images are useful for, among other things, detection of ischemia, 
cytotoxic edema, and demyehnization in the brain. Diffusion weighted images of the brain are 
complicated because the brain pulsates and to observe diffusion a series of images of the brain 
taken at different times are compared. Generally, very fast imaging techniques such as EPI are 
5 used in diffusion studies of the brain. Data acquisition for such imaging, with updating of a 
primary transform using a diffusion sensitive EPI sequence, in accordance with a preferred 
embodiment of the present invention can be useful in improving the quality of diffusion 
sensitive MRI imaging of the brain and other pulsating and moving organs. 

There is therefore provided in accordance with a preferred embodiment of the present 
10 invention a method for producing images of an object by MRI imaging comprising: acquiring a 
first set of data for the object at a first raster of coordinate points, in a first region in a k-space, 
using a first MRI imaging sequence; acquiring a second set of data for the object at a second 
raster of coordinate points in a second region in k-space using a second MRI imaging sequence, 
wherein the second region in k-space is included in the extent of the first region in k-space; 
15 using data from the second set of data to replace data from the first set of data thereby to 
produce a third set of data; and constructing a real space image of the object fi-om the third set 
of data. 

Preferably the real space image has substantially the same spatial resolution as a real 
space image constructed fi-om the first set of data. The second region in k-space preferably 
20 includes the origin of k-space. Preferably, the first region in k-space includes the origin of k- 
space. 

The second region in k-space is preferably significantly smaller than the first region in 
k-space. Preferably, the ratio of the size of the second region in k-space to the size of the first 
region in k-space is substantially equal to R^/(L-AL)^, where R is the resolution of an image of 
25 the object constructed from the first set of data, is substantially equal to the volume in real 
space of a feature observ^ed in the object and a change in the volume of the feature observed is 
substantially equal to [L^ - (L-AL P], and wherein n is equal to 2 and 3 for two dimensional 
and three dimensional imaging respectively. 

In some preferred embodiments of the invention data in the third set of data is smoothed 
30 to remove discontinuities. 

Preferably the time it takes to acquire the second set of data is substantially less than the 
time it takes to acquire the first set of data. 

10 
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The magnetic field strengths used in the first and second MRI imaging sequences are 
preferably within established medical safety guidelines. Preferably, the time derivatives of the 
magnetic fields used in the first and second MRI imaging sequences are within established 
medical safety guidelines. The second MRI imaging sequence preferably, uses substantially the 
5 same magnetic gradient field strengths as the first MRI imaging method. Preferably, the time 
derivatives of the magnetic fields used in the first and second MRI imaging sequences are 
substantially the same. Preferably, the second method uses substantially the same receiver 
bandwidth as the first MRI imaging sequence. 

In some preferred embodiments of the invention the first MRI imaging sequence is a 
10 gradient echo (GE) sequence or a spin echo (SE) sequence. The second MRI imaging sequence, 
in some preferred embodiments of the invention, is an EPI MRI imaging sequence. 

In some preferred embodiments of the invention the second MRI imaging sequence is 
T2 sensitive, in some others the second MRI imaging sequence is T2* sensitive. 

In some preferred embodiments of the present invention the imaging is diffusion 
1 5 weighted. 

Some MRI imaging methods in accordance with the present invention are sensitive to 
motion of parts of the body such as cardiac contractiUty, peristalsis, and respiratory motion. . 

In some preferred embodiments of the invention the MRI imaging method is sensitive 
to physiological changes such as those resulting fi-om neural activity, blood flow or the 
20 absorption of substance by different tissues in the body. 

In some preferred embodiments of the present invention, data in the first set of data that 
is not replaced is data that has remained substantially unchanged during the elapsed time 
between the acquisition of the first set of data and the acquisition of the second set of data and 
data that is replaced is data that has substantially changed during the elapsed time. 
25 Alternatively, or additionally, data in the first set of data that is replaced is data that defines 
larger features in the object and data that is not replaced is data that defines details of these 
features. 

BRIEF DESCRIPTION OF FIGURES 

The invention will be more clearly understood by reference to the following description 
30 of preferred embodiments thereof in conjunction with the figures, wherein common elements 
which appear in the figures are labeled with the same numeral, and in which: 

Figs. lA and IB show respectively, a conventional planar echo imaging sequence (EPI) 
and a raster of points in k-space at which the sequence samples data for a k-transform of an 
object; 
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Fig. 2 is a schematic of a hybrid two dimensional MRI imaging technique, in 
accordance with a preferred embodiment of the present invention, for the rapid production of a 
sequence of real time two dimensional images of an object; and 

Figs. 3 A-3D schematically illustrate elements of a hybrid imaging process for the rapid 
5 production of a sequence of real time three dimensional images of an object in accordance with 
a preferred embodiment of the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
A simple formula can be used to relate important parameters of MRI imaging sequences 
and clarify their significance to preferred embodiments of the present invention. 

10 The relationship KR~ 1 used above, can be written, yGT^R = 1, since K = yGT^, 

where y is the gyromagnetic ratio of the imaging nucleus, G the magnetic read gradient, and 
T/^ the data acquisition time in the MRI imaging sequence. Resolution R may also be written 
as the field of view of the image, FOV, divided by the number of pixels in the MRI image. 
Then KR = yGTa FOV/N = yGTo FOV - yG(l/(2BW)) FOV = 1, where Tq is the dwell time 

1 5 per point in k-space at which S( k ) is evaluated (Tq = T^/N, since the number of pixels in the 
real space image equals the number of data points taken in k-space), and BW is the bandwidth 
of the receiver (=1/(2Te))). 

Consider a conventional GE or SE sequence and assume that the imaging nucleus is the 
water proton, that a slice having an area 20cm x 20cm is imaged, and that the image has 256 x 

20 256 pixels {Le. R = 0.8 mm). Then y = 4257 Hz/gauss, FOV = 20, N - 256, and CTq ^ 1 1.75 
X 10-6. Assuming that a typical read time (i.e. T/^) for the GE or SE sequences is 3.072msec 
for a read line in k-space with 256 points, then Tq = 12 )j,sec, and G = 0.98 gauss/cm == 9.8 
mT/m. The net dwell time for all 256 x 256 points in k-space at which data is taken for the slice 
is about 1 second. However, the total image acquisition time, Tj, is equal to 256Tr, (for the 

25 case where each repetition of the sequence accumulates one read line of data) which is 
considerably longer than the total dwell time if Tr is on the order of Tj . For brain white matter 
for example, Tj is about 0.5 sec, which would lead to a total image time Tj, for a single slice of 
brain tissue, using conventional GE and SE sequences, of approximately two minutes. The 
resolution of the GE and SE sequences might be acceptable for most real time imaging but the 

30 imaging time is too long. Not only is the total image time Tj much too long for rapid imaging 
of important physiological processes, but even the net dwell time is too long ( e.g. it's too long 
to image a heart beat). 

12 
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In comparison, consider an EPI sequence that acquires data at all 256 x256 data points 
in k-space for the slice in a single sequence, and that the imaging time is an acceptable, 100 
msec. For this case the dwell time becomes 1.5|isec, the gradient G becomes 78.4mT/m, and 
the time rate of change of the read gradient field is on the order of 64 times that used in the GE 
5 and SE sequences discussed in the previous example. This time rate of change of field is large 
aud potentially harmful to biological tissue. 

In addition, whereas in the case of the SE and GE sequences receiver bandwidth is 
= 42kHz, for the EPI sequence it is 336kHz. Therefore, although the EPI sequence can perform 
fast enough, it has serious technical handicaps, and as noted above, it does not generally have 
1 0 the resolution required. 

Now consider using fast EPI to update only the central 64 x 64 k-space data points for 
the slice, in accordance with a preferred embodiment of the present invention. Assume a total 
net dwell time of 50msec. Then Tj) is again 12|j,sec and bandwidth is = 42kHz. Using the 
relationship yGTj) FOV = 1, G has the value 9.8 mT/m, the same value used by the GE and 
15 SE sequences. Rise and fall times of the gradient field can be an acceptable 100 to 150(asec for 
read field gradient blips in the EPI sequence (see Fig. 2A). Such rise and fall times extend the 
acquisition time of the 64x64 data points from the net dwell time of 50msec to between 60 and 
70msec. 

It is therefore seen that a hybrid imaging process in accordance with a preferred 
20 embodiment of the present invention, that preferably uses a GE or SE imaging sequence to 
acquire a high resolution primary k-transform of an object imaged and a fast EPI imaging 
sequence to update part of the primary k-transform can produce real space images of the object 
with substantially the rapidity of the fast EPI sequence and substantially the resolution of the 
high resolution GE or SE sequence. It is also seen that the hybrid imaging process in 
25 accordance with the present invention can do this with conventional gradient strengths and 
receiver bandwidths, Le. with conventional GE or SE equipment. 

The imaging sequences used in a hybrid imaging sequence in accordance with a 
preferred embodiment of the present invention can be tailored to a type of image and type of 
image contrast desired. For example, for diffiision imaging of the brain a diffusion pulse train 
30 can be added to the imaging sequence and preferably updated with snapshot EPI techniques. In 
some preferred embodiments of the invention the sequence provides T2 contrast. In other 
preferred embodiments the sequence provides T2* contrast. Other preferred embodiments of 
the present invention can be made sensitive to temporal changes in the body in order to image 

13 
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peristalsis or cardiac contractility, or physiological changes resulting from neural activity. Still 
other preferred embodiment can be tailored for use with contrast agents such as Gd-DTPA or 
for imaging changes in tissues due to oxygen uptake. 

Fig. 2 shows, schematically, a two dimensional MRI hybrid imaging process 60, in 
5 accordance with a preferred embodiment of the present invention, for rapid sequential updating 
of MRI images in a real time study of an object. 

A primary two dimensional k-transform, So(kx,ky), of the object is acquired at a time 
tQ, over a primary raster 62, of points 64, in k-space. Preferably, primary raster 62 comprises 
256 equally spaced read lines 66, having 256 equally spaced read points each. Preferably, the 
10 perimeter of primary raster 62 is a square 67, whose center is at the origin of k-space. 
So(kx,ky) is calibrated and phase corrected using methods known in the art so that the k-space 
origin is correctly centered. A primary image IMq of the object is constructed by Fourier 
transforming primary k-transform So(kx,ky). 

Primary k-transform So(kx,ky) is preferably acquired using a high resolution MRI 
15 imaging sequence such as preferably, a GE or SE sequence of the type shown in Figs. lA and 
IC respectively. 

Phase correction of the primary k-transform should not be required if a SE sequence is 
used. However if a GE technique is used, phase correction can preferably be carried by 
acquiring calibration data for So(kx,ky) at read points along a single kx read line without phase 

20 encoding along ky. The calibration data is Fourier transformed and the phase of each pixel 
along the x axis calculated. A complete set of data for So(kx,ky) is then acquired and Fourier 
transformed for the kx variable only, to get an interim function IF(x,ky). Interim function 
IF(x,ky) is a function of the real space variable x and the k-space variable ky. The phases 
calculated from the calibration read line are complex subtracted for each pixel, point by point, 

25 for every read line of the function IF(x,ky). and the phase corrected IF(x,ky) is Fourier 
transformed for the ky variable to get a phase corrected primary image IMq. Primary image 
IMq is then Fourier transformed back into k-space to yield a primary k-transform. 

At a subsequent time, tj, an update k-transform, Su(kx,ky), of the object is acquired at 
an update raster 68, of read points 70. Preferably, update raster 68 comprises 64 read lines 72, 

30 having 64 equally spaced read points each. Preferably, the spacing between read Hnes 72 and 
read points 70 in update raster 68 are the same as the spacing between read lines 66 and read 
points 64 in primary raster 62. Preferably, the perimeter 69, of update raster 68 is a square 
whose center is at the origin of k-space. The area of update raster 68 is therefore congruent 
with a central square area 74, of primary raster 62. which is 64 read lines high by 64 read 

14 
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points wide and read points in raster 68 are congruent with read points 64 in central area 74. 
Like So(kx,ky), update k-transform Su(kx,ky) is calibrated and phase corrected using methods 
known in the art so that the k-space origin is correctly centered. 

Values of So(kx,ky), at read points 64 in central area 74 of primar\' raster 62 are then 
5 replaced with values calculated from update image Su(kx,ky) to produce a new k-transform of 
the object. The new k-transform is Fourier transformed to produce an updated image, IMl, of 
the object at time tl. The process is repeated again at a time t2, and subsequently as many 
times as needed and conditions allow. The time between images can be as short as the update 
time, which is on the order of from 20- 100msec, if not limited by other parameters such as Tj 
10 recovery time. 

Prior to replacement, amplitudes of Su(kx,ky) are calibrated to amplitudes of So(kx,ky). 
If there is a mismatch between the spacing of read lines and read points in update raster 68 and 
central area 74 of raster 62 then measured values of Su(kx,ky) are used to interpolate 
appropriate values of Su(kx,ky) to replace values of So(kx,ky). 

15 Preferably, read points 70 along the inside edge of the area of update raster 68 fall 

within an area of primary raster 62 where the k-transform of the object does not change 
significantly during the real time study of the object. This facilitates the calibration of 
amplitudes between Su(kx,ky) and So(kx,ky) and the interpolation of values of Su(kx,ky) 
when necessary in the case where read lines between primary raster 62 and update raster 68 get 

20 misaligned. Aspects of matching data and smoothing of transitions at the boundaries between 
two sets of related k-transform data are discussed in a patent application titled, "Removing 
Discontinuities in MRI K-Space Data", filed in the US Patent Office by a same applicant as an 
apphcant of this application, on September 14, 1998, the disclosure of which is incorporated 
herein by reference. 

25 It should be realized that the full primary image can be updated as needed or desired. 

An indication as to when this is required can be inferred from the quahty of the images 
constructed from the updated k-transforms. 

Three dimensional hybrid imaging processes for rapid sequential imaging, in 
accordance with a preferred embodiment of the present invention, are based on the same 
30 principles as the two dimensional processes described above. 

Figs. 3A-3D illustrate elements of a three dimensional hybrid imaging process for 
imaging a thick slice of an object in accordance with a preferred embodiment of the present 
invention. 
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A high resolution three dimensional primary k-transform So(kx,ky,k2) of the object is 
acquired at read points in a primary three dimensional raster in k-space. Preferably, the 
acquisition is made using a three dimensional GE imaging sequence 80 of the type shown in 
Fig. 3A. In imaging sequence 80 a broad band selection RF pulse 82 and selection gradient 84 
5 select a thick slice of the object for imaging. Phase encoding is done in both the y and z 
directions using y phase gradient 88 and z phase gradient 90 respectively. Y phase gradient 88 
is step increased by a fixed y phase increment 89 for each read line scanned and z phase 
gradient 90 is step increased by a fixed z phase increment 91 after all the lines in a read plane 
are scanned. Step increases for the y and z phases are performed between repetitions of GE 

10 sequence 80. Data is acquired during x read gradient 92 by sampling data echo 96. Gradient 94 
which precedes read gradient 92initializes the starting k^ value for each read line and centers 
data echo pulse 96 at the center of read gradient 92. 

Preferably, the raster of read points at which So(kx,ky,k2;) is evaluated comprises L x N 
x M read points in k-space. Preferably the read points are distributed on L identical parallel 

1 5 plane surfaces equally spaced along the z axis, each plane surface having preferably the same 
number N x M read points. Preferably, the N x M read points in a plane lie on N read lines of 
M read points each, parallel to the kx axis. Preferably, the distance between adjacent read lines 
on a plane is the same for all read lines. The distance between a read point and an adjacent read 
point on a same read line is preferably the same for any read point in the raster. 

20 Reference is now made to Fig. 3B. Preferably, the L plane surfaces of the primary raster 

(not shown) are contained in a primary k- volume 100 which is rectangular and centered on the 
k-space origin. 

Following the acquisition of the primary k-transform, three dimensional update k- 
transforms, Su(kx,ky,k2), are preferably evaluated at a reduced number of read points (not 
25 shown) in a rectangular update k-volume 102 contained inside primary k-volume 100. Update 
k- volume 102 is preferably centered on the origin of k-space and much smaller than primary k- 
volume 100. 

As in the two dimensional case, data from the three dimensional update k-transform 
Su(kx,ky,k2) replaces data in the three dimensional primary k-transform So(kx,ky,k2) in order 
30 to sequentially image the observed object in real time. The relative size of the update k-volume 
to the primary k-volume (which is generally the ratio of the number of data points in the update 
k-transform to the number of data points in the primary k-transform) is preferably, substantially 
equal to r3/(L-AL)3 where R is the resolution of the real space images L is the characteristic 
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size of the features in the object observed and AL is the characteristic size of the changes in 
these features. 

Again, as in the two dimensional case, the reduced number of data points required for a 
three dimensional update makes the use of very fast three dimensional MRI imaging sequences, 
5 such as three dimensional EPI, advantageous. As a result, updates can be made very rapidly 
and MRI imaging procedures, in accordance with a preferred embodiment of the present 
invention, can provide rapid three dimensional real time sequential imaging with high 
resolution. 

A three dimensional EPI sequence 136, preferably used for accumulating data for 

10 update k-transforms is shown in Fig. 3C. For clarity and ease of exposition, EPI sequence 120 
refers to an update volume enclosing a raster of 4x4x4 read points, Le. 4 read planes (L=4) 
having 4 read lines (N=4) with each read line having 4 read points (M=4). 

Following a selection pulse 122, with a flip angle a, a read gradient blip train 124 
comprising 16 read blips 126 is appHed to the x axis. A phase gradient blip train 128 

15 comprising 12 phase blips 130 is applied to the y axis and a phase gradient blip train 132 
comprising 3 phase blips 134 is applied to the z axis. Blip trains 124, 128 and 132, begin with 
negative initialization gradients 136, 138 and 140 respectively, which initialize the coordinates 
of the starting read point of the data scan. 

Since there are 1 6 read lines in the raster, and all the read points on a read line are read 

20 during one read gradient blip 126, there are 16 read gradient blips 126 in read gradient blip 
train 124. Each read blip is lasts for a time Ty^, All except every fourth read gradient blip 126 
in blip train 124 is followed with a y phase gradient blip 130 which is centered at the zero 
crossing point between read blips 126. Each y phase gradient blip 130 lasts for a time Ty, and 
increases ky by yGyTy, thereby advancing the k-scan to the next read line in the read plane. 

25 Every fourth read gradient blip 126 is followed by a z phase gradient blip 134 which lasts for a 
time T2 and increases kz by 70^X2 to advance the k-scan to the next read plane. Z phase 
gradient blips 134, when they occur, Uke y phase gradient blips 130, are centered on the zero 
crossover points between read blips 126. Y phase gradient blips 130 change sign following 
each z phase gradient blip 134 because from one read plane to the next, the direction of the k- 

30 scan reverses in the y direction. Data is recorded from a train 140, of echo pulses comprising 
data echo pulses 142. There is an echo pulse for each read gradient blip 126. The maximum 
amplitude echo pulse generally occurs at the origin of k-space unless the echoes have been 
significantly reduced by transverse T2 or T2* dephasing, 
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Total data acquisition time Tj^ for an update volume using an EPI scan of this type is 
equal to NxLxTfyj- With an amplitude of about 9.8 mT/m and rise time of 116|isec for read 
gradient blips 139, then for a read line having 16 read points and dwell time of 12|asec per read 
point Tf^ = 424|j.sec- Data for a complete update raster of 16x16x16 read points is then 
5 acquired in 256Tjvi= 108msec. 

Fig. 3D shows a k-space trajectory 154, traversed by 3 dimensional EPI sequence 120 
as it samples data at 4x4x4 read points 156, in an update raster 154. Raster 154 comprises 4 
read planes 158 separated by a distance in k-space equal to 762X2, where y is the gyromagnetic 
ratio of the imaging nucleus, G2 is the amplitude of z phase blips 143 and T2; is the duration of 

10 phase blips 143. The relative displacement of the read planes along the z axis is noted by the 
numeral at the lower right of the read plane. Read plane 158 labeled with "0" coincides with 
the kxky plane, read plane 158 labeled "-1" is displaced along the negative k2 axis by YG2T2. 
Each read plane 158 has 4 read lines 160 of length yGxT^ separated by yGyTy, where G^, 
Tm, Gy, and, and Ty are the gradient blip amplitudes and times of read gradient blips 139 and 

15 phase blips 141 respectively. The data scan starts at read point 162 and proceeds along 
trajectory 152 in the direction of arrows 164. 

Similarly to the two dimensional case, for three dimensional rapid sequential imaging 
of an object, in accordance with a preferred embodiment of the present invention, spacing 
between read planes, read lines and read points in an update raster are preferably identical to 

20 spacing between read planes, read lines, and read points in the primary raster for which the 
update is intended. Preferably, the update raster includes read points in k-space for which the 
values of a primary k-transform of the object do not change much. This facilitates calibrating 
update k-transform data to the primary k-transform data and smoothing transitions between 
update data and data retained from the primary k-transform that is used for imaging at times 

25 subsequent to the time when the primary k-transform is acquired. 

It should be realized that update rasters do not necessarily have to occupy the same 
volume of k-space for each update. As changes in an object imaged accumulate over time, data 
that was initially sufficiently accurate to be retained from a primary k-transform of the object 
and used for imaging at one time, might degrade at a later time to such an extent that it should 

30 be replaced. In such cases the volume of a raster of read points might be enlarged or moved or 
otherwise altered to satisfy and anticipate needs. 

The present invention has been described using a non limiting detailed description of a 
preferred embodiment thereof. Variations of the embodiment described will occur to persons 
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of the art. The detailed description is provided by way of example and is not meant to limit the 
scope of the invention which is Umited only by the following claims: 
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CLAIMS 

1 . A method for producing images of an object by MRI imaging comprising: 

acquiring a first set of data for the object at a first raster of coordinate points, in a first 
region in a k-space, using a first MRI imaging sequence; 
5 acquiring a second set of data for the object at a second raster of coordinate points in a 

second region in k-space using a second MRI imaging sequence, wherein the second region in 
k-space is included in the extent of the first region in k-space; 

using data firom the second set of data to replace data firom the first set of data thereby 
to produce a third set of data; and 
10 constructing a real space image of the object fi-om the third set of data. 

2. A method according to claim 1 wherein the real space image has substantially the same 
spatial resolution as a real space image constructed firom the first set of data. 

15 3. A method according to claim lor claim 2 wherein the second region in k-space includes 
the origin of k-space. 

4. A method according to any of the preceding claims wherein the first region in 
k-space includes the origin of k-space. 

20 

5. A method according to any of the preceding claims wherein the second region in k- 
space is significantly smaller than the first region in k-space. 

6. A method according to claim 2 wherein the ratio of the size of the second region in k- 
25 space to the size of the first region in k-space is substantially equal to RW(L-AL)^, where R is 

the resolution of an image of the object constructed from the first set of data, is substantially 
equal to the volume in real space of a feature observed in the object and a change in the volume 
of the feature observed is substantially equal to [L^ - (L-AL )^], and wherein n is equal to 2 
and 3 for two dimensional and three dimensional imaging respectively. 

30 

7. A method according to any of the preceding claims wherein data in the third set of data 
is smoothed to remove discontinuities. 
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8. A method according to any of the preceding claims wherein the time it takes to acquire 
the second set of data is substantially less than the time it takes to acquire the first set of data. 

9. A method according to any of the preceding claims wherein the magnetic field strengths 
used in the first and second MRI imaging sequences are within estabhshed medical safety 
guidelines. 

10. A method according to any of the preceding claims wherein the time derivatives of the 
magnetic fields used in the first and second MRI imaging sequences are within estabhshed 
medical safety guidelines. 

11. A method according to any of the preceding claims wherein the second MRI imaging 
sequence uses substantially the same magnetic gradient field strengths as the first MRI imaging 
method. 

12. A method according to any of the preceding claims wherein the time derivatives of the 
magnetic fields used in the first and second MRI imaging sequences are substantially the same. 

13. A method according to any of the preceding claims wherein the second method uses 
substantially the same receiver bandwidth as the first MRI imaging sequence. 

14. A method according to any of the preceding claims wherein the first MRI imaging 
sequence is a gradient echo (GE) sequence. 

15. A method according to any of claims 1-13 wherein the first MRI imaging sequence is a 
spin echo (SE) sequence. 

16. A method according to any of the preceding claims wherein the second MRI imaging 
sequence is an EPI MRI imaging sequence. 

17. A method according to any of the preceding claims wherein the second MRI imaging 
sequence is T2 sensitive. 
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18. A method according to any of the preceding claims wherein the second MRI imaging 
sequence is T2* sensitive. 

19. A method according to any of the preceding claims wherein the imaging is diffusion 
weighted. 

20. A method according to any of the preceding claims wherein the imaging is sensitive to 
motion of parts of the body such as cardiac contractility, peristalsis, and respiratory motion. 

21. A method according to any of the preceding claims wherein the imaging is sensitive to 
physiological changes such as those resulting jfrom neural activity, blood flow or the absorption 
of substances by different tissues in the body. 

22. A method according to any of the preceding claims wherein data in the first set of data 
that is not replaced is data that has remained substantially unchanged during the elapsed time 
between the acquisition of the first set of data and the acquisition of the second set of data and 
data that is replaced is data that has substantially changed during the elapsed time. 

23. A method according to any of the preceding claims wherein data in the first set of data 
that is replaced is data that defines larger features in the object and data that is not replaced is 
data that defines details of these features. 
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